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PACS 61.20.Qg - Structure of associated liquids: electrolytes, molten salts, etc. 
PACS 82 . 45 . G j - Electrolytes 
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Abstract - A Poisson-Boltzmann approach is used to determine the double-layer integral and 
differential capacitances in a finite-length situation for an electrolytic cell. By means of simple 
analytical calculations, it is shown how these quantities exhibit the bell-like and the camel-like 
shapes as a function of the applied voltage, when the thickness of the sample or the surface charge 
are varied. The problem is formulated treating the bulk liquid as a member of a grand canonical 
ensemble in contact with blocking electrodes. As a consequence, an applied voltage dependent 
Debye's screening length arises as a fundamental length governing the electrical behavior of the 
system. 



The capacitance curves in ionic liquids have typically 
one (bell-likc shape) or two (canicl-likc shape) maxima, 
when plotted as a function of the applied voltage. In re- 
cent papers, this behavior has been explained as resulting 
from an interplay between the excluded volume effect and 
the effective electric constant [1] or as being essentially 
governed by a single parameter that represents the ratio 
of the bulk density of ions in the liquid to the maximal 
density in the double layer In a recent article, Hatlo 
et al. [3] considered the excess ion polarizability into the 
Poisson-Boltzmann theory in order to show that the de- 
crease in differential capacitance with voltage can be un- 
derstood in terms of thickening of the double layer due 
to ion induced polarizability holes in water. In [4] was 
proposed a theory for concentrated electrolyte solutions 
in the context of lattice model and Fermi statistics to con- 
sider ions size effect for the diffuse double layer. It seems 
that all these important approaches invoke, in some de- 
gree, the necessity to go beyond the Poisson-Boltzmann 
framework to face the very rich behavior of the experi- 
mental data regarding the double-layer capacitance in a 
typical electrolytic cell. In addition to what is affirmed in 
these papers, here we show that by varying the thickness 
of the sample, some essential features of these results, like 
the camel- and bell-shaped behavior of the capacitance of 
an electrolytic cell as a function of the applied voltage, are 



obtained. The model permits to show, in a simple man- 
ner, how the capacitance (integral or differential) curve of 
a finite sample of an insulating medium containing mo- 
bile ions can be easily determined in the framework of 
the Poisson-Boltzmann theory even for point-like parti- 
cles. Moreover, we show also how the thickness of the 
sample is responsible for the shift in the maxima of the 
capacitance that can be found at low voltages. Recently, 
in Ref. [5] the nonlinear Poisson-Boltzmann theory was 
invoked to investigate the behavior of a non-conductive 
quasi-planar lipid membrane in an electrolyte and in a 
static (dc) electric field. 

We consider an insulating medium containing ions, with 
uniform dielectric coefhcient e, in the shape of a slab, with 
two uniform flat surfaces in the x — y plane, separated by 
a distance d, located at z = ±(i/2 where blocking elec- 
trodes are placed. It is supposed, for simplicity, that all 
the physical quantities entering in the calculation are only 
z— dependent. The fundamental equations have been pre- 
sented in [6ll7] but they will be summarized here for clarity. 
In the absence of external field the medium is global and 
locally neutral such that the densities of positive and neg- 
ative ions are equal to the initial density, i.e. n±{z) = hq. 
In the presence of a difference of potential U (supplied by 
an external agent) such that V{zkd/2) = ztU/2 the equilib- 
rium distribution of charges is governed by classical statis- 
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tics, i.e., n±{z) = ne^^'^^\ in which = qy{z)/kBT 
is the electric potential in ksT/q units, where ks is the 
Boltzmann constant, T is the absolute temperature, and 
q is the magnitude of the elementary charge. Notice that, 
in the present formalism n is the density of ions in the 
liquid phase, to be determined by imposing the conserva- 
tion of the number of ions. This approach is equivalent 
to working with a fluctuating number of particles in the 
bulk and an unknown chemical potential that is formally 
determined by the conservation of the number of particles. 
It is just this condition that is responsible for the incor- 
poration of a voltage dependence to the Debye's screening 
length, as shown below. The spatial dependence of the 
electrical potential is governed by the Poisson's equation, 
which, after recognizing that the density of electric charge 
is given by p{z) = q[n+{z) — n-{z)\, can be put in the 
form 
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is a length depending on the applied voltage and connected 
with the Debye's screening length Aq = eksT / {2q^nQ). 
The nonlinear eq. ([T]) has to be solved with the bound- 
ary conditions '0(±(i/2) = ±u = qU/{2kBT), and can be 
integrated to give 



y^cosiri/r+T, 
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where k is an integration constant. The conservation of 
the number of particles is given by 



nod= / n+{z)dz^—r=n 
and can be rewritten as 
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whereas eq. ([2]) becomes 
1 
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which, together, imply that 
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The fundamental equations have been established and 
the procedure is now straightforward. Indeed, I{k, u) and 
J(fc, u) can be expressed in terms of elliptic integrals of 
first and second kind and, for a given applied voltage [/, 
from eq. ([6]), one obtains k{U) for a sample of definite 
thickness d and a Debye's screening length Ad- Once this 
is done, the effective (voltage dependent) Dcbye screening 
length can determined by the expression ([5]) . Finally, with 
these quantities determined as a function of U, the elec- 
tric field profile is then easily obtained from eq. ([2]), i.e., 
E{z) = -{ksT/q) d^j/dz. 
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Fig. 1: Integral double-layer capacitance vs. the bias voltage, 
U, for different thicknesses of the sample: d = 10 fim (thick 
solid), d — 8 fim (dashed-dotted) , d = 5 fim (dashed), d = 3 
(dotted), and d = 1 fj,m (thin solid). The curves have been 
drawn for Ad = 0.5 /xm. 



To obtain simple expressions for the double-layer capac- 
itance, it is necessary to calculate the net charge on the 
surfaces, given by cr = eE(z ~ — d/2). In this manner, the 
difi^erential double-layer capacitance is given by 



I{k, u)\/cos\iu -I- k 



da C ^ 
dU du 

in which Cq = e/2d and © and ([S]) have been used. In the 
same manner, the integral double-layer capacitance will be 
given by 
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In fig. ([TJ , the integral capacitance is shown as a function 
of the external voltage. The camel-like shape is evident, 
but it is strongly dependent on the ratio between the thick- 
ness of the sample and the Dcbye length Ad ■ Indeed, when 
~ Ad, the maximum is found for t/ = (see, e.g., the 
thin solid line). The camel-like shape is found if c? > Ad. 
Likewise, the differential double-layer capacitance, shown 
in fig. [5J exhibits similar behavior with the thickness. In 
fig. [21 the integral double-layer capacitance is shown as a 
function of the charge a for different thicknesses of the 
sample. The analysis presented above is based on the 
conservation of the number of particles in the capacitor, 
eq. (3). This means that, at first sight, we are considering 
a closed system, where the number of ions is fixed by its 
initial volume, as in a sealed condenser. In this manner, 
when the system is subject to higher voltages the amount 
of net charge accumulated in each double layer will even- 
tually stop growing with further increase of voltage, and 
the capacitance will decrease. For this reasons we observe 
that in larger systems the falling branches of capacitance 
emerge at larger voltages. The application of our model 
to real experiments, where the electrodes are embedded 
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Fig. 2: Differential double-layer capacitance vs. the bias volt- 
age, U, for different thicknesses of the sample according to the 
convention adopted in fig. [1] 



into a sea of electrolyte, and can suck ions from the bulk, 
requires some caution. In fact, in this case the system is 
no longer closed, and it can exchange particles with the 
reservoir. However, we are working here in the slab ap- 
proximation, i.e., the electrolyte is limited by two flat, 
infinitely large, electrodes. This means that our sample is 
isolated, as the semi-infinite samples considered in various 
recent works [HlHlin] • It follows that our theoretical predic- 
tions can be compared with the ones reported in |2j[8n9j. 
Since the quoted theoretical model [2l[8l|9] have been used 
to interpret experimental data relevant to an open con- 
denser, where the electrodes can adsorb ions coming from 
outside of the sample, we can conclude as follows. In real 
case, the systems are not isolated, but the dominant role 
on the electrical response of the electrolytic cell is played 
by the ions contained in the volume limited by the elec- 
trodes. In the opposite case, the theoretical description 
proposed by us, as the one following from the models dis- 
cussed in [21IH1IS], does not work, and the description has 
to include also the lateral ionic current. However, at the 
present, due to the good agreement between theoretical 
predictions and experimental data it seems that the effect 
related to these horizontal current can be neglected. 

The behaviors shown in fig. 1, fig. 2 and in fig. 3 are 
consistent with the results found in [S] as a function of the 
concentration and with the Monte Carlo simulation in [5] . 
The simulation were performed in the grand canonical en- 
semble to remove uncertainties around the determination 
of the bulk concentration. Likewise, in our system, for a 
fixed value of no, the actual bulk concentration depends 
on n which, in turn, is fixed by imposing the conservation 
of the number of particles. This means that our approach 
treats the liquid system as a member of a grand canon- 
ical ensemble as well, i.e., as a system having a variable 
number of particles and, from the point of view of the 
ensemble theory, it is equivalent to have as macroscopic 
variables T, V, /i, where n is the chemical potential. The 
complete system, i.e., the liquid medium plus the walls 
(solid phase) is, instead, characterized by the equilibrium 



Fig. 3: Integral double-layer capacitance vs. surface charge, 
for different thicknesses of the sample as in figs. [1] and (2] 



variables T, V, N, i.e., it can be treated in the framework 
of the canonical ensemble. In addition, we notice that a 
bell-like shape curve similar to the one shown in figs. ([T|) 
and ([2]) has been found in the framework of a mean-field 
approach at higher concentration |10j and by means a 
Landau-Ginzburg-type continuum theory which takes 
into account the exclusion volume corrections. At low 
concentrations, a camel-like shape capacitance curve was 
reported in [l][2]- In our approach, illustraded here for 
parameters more appropriated to ordinary electrolytes, 
these behaviors are found when the ratio d/Ao varies. In 
fact, when d/Xo > 1 (which corresponds to a more di- 
luted system), we obtain the camel-shape behavior. For 
d ~ Xdi the model predicts a bell-shape behavior. Finally, 
in fig. m the integral double-layer capacitance as a func- 
tion of the external voltage is shown in the limiting situa- 
tion for which A — > Ad , and exhibits the expected classical 
Gouy-Chapman behavior [T2UI5] . In this case, the voltage 
dependence introduced in X{U) by means of the condition 
of conservation of number of particles is lost, as can be 
deduced by comparing the behaviors shown in figs. [T] and 
m In fig.m the effective Debye's length depends on U and 
we can found both, the bell- and camel-shape behavior as 
d/Xo increases. In fig. SI A = Xd, and the curve presents 
the Gouy-Chapman behavior, i.e., it is minimum at the 
point of zero charge. 

In the high voltage limit, i.e., it 1, it is possible 
to show that the electric potential across the sample be- 
comes [5]: 



ip{z) = 2u- 



2 sinh — 

u \ d 



where D = d/(2AD)- The differential double-layer capaci- 
tance is now given by 
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Cd = 2Co 

whereas the integral capacitance becomes 
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Fig. 4: Integral double-layer capacitance vs. the applied volt- 
age, for d = IQ fim and Ad = 0.5 ^m, when A = Ad. 
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In the limiting situation of small applied voltage, i.e., u <C 
1, the electric potential may be approximated as [6]: 



sinh(z/AD) 



sinh(d/2AD)' 

and the net charge on the surface is given by cr = 
(e/2AD)C^cothD, and, in this case, both the differential 
and the integral capacitance coincides, namely 

Cd^a = CGc cothD^Ccc, 

in which Cqc = (£/2Ad) is the linear Gouy-Chapman or 
"Debye" capacitance [TB], as expected because, typically, 
I? » 1, and cothZ? w 1. 

Actually, in the limit of very large separation between 
the electrodes, the sample can be considered as formed by 
two half spaces, as shown in [6], and the electric potential 
may be written as ip{z) = ip+^z) + where 



?A± = 2 In 



1 ± 76 



{zTd/2}/\D 
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with 7 = tanh — 
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By using eqs. ([7]) and ([5]), this exact limiting expression 
leads to 



2^*00 cosh 



and Ci = — 
2 



sinh(l) 
^Cgc- 
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for D ^ 1. Thus, for very large separation between the 
electrodes, one concludes that the capacitances are inde- 
pendent of this separation. Moreover, for it = 0, both the 
differential as well as the integral capacitances reduce to 
the half of the Couy-Chapmann capacitance correspond- 
ing to a semi-infinite sample. This is not surprising, be- 
cause for very large sample one should not expect to ob- 
serve a difference between the grand-canonical and canon- 
ical ensembles, since most of the ions are coming from the 



Fig. 5: Differential double-layer capacitance vs. the applied 
voltage, for d = 10 and Ad = 0.5 ^m, for three different 
temperatures: T = 250K (dashed), T = 290K (dotted) and 
T = 298.15 K (solid). 



bulk. To complete the analysis, in fig. ([5]), the differen- 
tial double-layer capacitance is shown as a function of the 
applied voltage for different values of the temperature. It 
diminishes with increasing temperature as found in other 
approaches [9lfT6]. 

Even if it is true that electrochemical process, such as 
reactions, are determined by the ion profiles [17], and the 
capacitance is some kind of integral of ion profiles, its de- 
termination is crucial for investigating the response of the 
electrolytic cell to an external applied voltage. The very 
rich behavior of the capacitance is a confirmation of its im- 
portant role in the electrical behavior of the cell. In this 
sense, the efforts concentrated in the understanding of this 
complex behavior are justified. Here, we have presented 
a simple analysis giving a complete profile for the behav- 
ior of the double-layer capacitance in a finite-length situ- 
ation in the framework of the Poisson-Boltzmann statis- 
tics. For explaining the transition from the bell-shape to 
the camel-shape behavior, it is crucial to describe in the 
correct manner how the bulk concentration of particles 
are varying in the presence of the applied voltage. In the 
formalism we have presented, this is done by considering 
a fluctuating bulk density which is fixed by the imposi- 
tion of the conservation of the number of particles. More 
precisely, even working with blocking electrodes, this for- 
malism permits one to treat the charges accumulated at 
the electrodes as pertaining to the surface of the system, 
instead of remaining in the bulk. As a consequence, in the 
system appears a Debye's screening length that depends 
on the bias voltage. Only in the absence of external field, 
this quantity coincides with the usual Debye's screening 
length. It is then possible to account for the effect of the 
external voltage by considering the role of two important 
lengths characterizing the system, namely, the thickness 
defining the geometry of the cell, and an "effective" elec- 
tric length strongly connected with the Debye's classical 
length. In this manner, for d ~ A^, the system presents 
a bell-shaped (a maximum at the point of zero charge) 
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behavior. As d increases, a shallow minimum appears at 
the point of zero charge and two other symmetric maxima 
appear in the region of low voltage (camel-shaped behav- 
ior). The shallow minimum evolves to a deep minimum 
when d ^ Xjj. As a final remark, it is worth mention- 
ing that all the more sophisticated approaches that have 
been proposed in the last years to face the problem have 
enlarged a lot the comprehension of the problem and new 
efforts have to be made to push further this comprehen- 
sion. However, a global understanding of the electrical 
behavior of the sample that takes into account the in- 
tegrated ion profiles can be easily achieved by means of 
a simple Poisson-Boltzmann approach formulated in the 
framework of the grand canonical ensemble. 
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